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Coulomb Effects in Nanoscale SINIS Junction 
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We study the system of two superconductors connected by a small normal grain. We consider the modifi- 
cation of the Josephson effect by the Coulomb interaction on the grain. Coherent charge transport through the 
junction is suppressed by Coulomb repulsion. An optional gate electrode may relax the charge blocking and 
enhance the current leading to the single Cooper pair transistor effect. Temperature dependences of critical 
current and of the minigap induced in the normal grain by the proximity to superconductor are studied. Both 
temperature and Coulomb interaction suppress critical current and minigap but their interplay may lead to 
the nonmonotonous and even reentrant temperature dependence. 

PACS: 73.21.-b, 74.50.+r, 74.45.+C 



Nanoscale SINIS junction consists of a small normal 
metallic grain connected to two superconductive leads 
by tunnel junctions. The Josephson effect in this sys- 
tem is provided by the Andreev reflection processes at 
the contacts: Cooper pairs enter the normal part of the 
junction and propagate as Cooperons. This allows the 
transport of Cooper pairs from one lead to another es- 
tablishing a supercurrent. Another manifestation of the 
Andreev mechanism is the appearance of a minigap in 
the spectrum of the normal grain: the proximity effect. 
This minigap is a result of a non-zero Cooper pairs den- 
sity come from superconductors. The mentioned effects 
rely on the phase conservation in the grain. Charging 
effects lead to fluctuations of the phase and break the 
coherent electron transport as well as the induced mini- 
gap. The interplay between proximity and charging ef- 
fects in a normal grain connected to one superconductor 
was recently studied in pQ . Here we apply the same for- 
malism to the system with two leads and consider the 
Josephson effect. We also study temperature depen- 
dence of the minigap and Josephson current, and find 
some unexpected reentrant behavior in a certain range 
of parameters. 

We consider tunnel junctions between normal grain 
and left (L) and right (R) superconductors characterized 
by large [in units of e 2 /h] normal-state conductances 
Gl,r 3> 1 and (geometric) capacitances C^r. The gate 
electrode is coupled to the grain by the capacitance C g . 
Mean level spacing in the grain 5 is the smallest energy 
scale of the system while the superconductive gap A 
in the leads and Thouless energy £rh of the grain are 
largest ones. We assume that Andreev conductances 
of both contacts are small, R < 1 (together with 



conditions Gz, r ^ 1 it means that our junctions con- 
tain many weakly-transparent channels, quantitatve es- 
timates will be provided below). The proximity and 
charging effects in the grain are characterized by the 
bare minigap width E g o = (Gz, + Gr)8/A and Coulomb 
energy E c = e 2 /(2C) with C = C L +C R + C g + AC 
being the total capacitance of the grain. Here AC = 
Ta(Gl + Gr) is the contribution to capacitance com- 
ing from virtual quasiparticle tunnelling . We assume 
that 6 « (Ego,E c ) < (A,E Th ). With the help of the 
dynamical (in imaginary time) sigma-model in replica 
space 0], and the adiabatic approximation for charging 
effects developed in PQ, we study the current-phase rela- 
tion of SINIS junction as well as the dependence of the 
critical current on temperature and gate voltage V g . 

Electronic properties of the normal grain are charac- 
terized by its Green function. To capture proximity in- 
duced correlations one uses the matrix Green function in 
Nambu-Gor'kov representation. The sigma-model oper- 
ates with the matrix field Q which apart from Nambu- 
Gor'kov structure carries two Matsubara energies and 
two replica indices. The standard Green function can 
be extracted from the diagonal in energies element of 
matrix Q by replica averaging. Charging effects are de- 
scribed by the fluctuating scalar field <f) corresponding 
to the electric potential and also carrying the replica 
index in the sigma-model formalism. The action for the 
two variables Q and <f> reads 



S[Q,cf>} =-- Tr [(et 3 + <f>)Q] 

l/T 

-^Tt[{G L Q L + G R Q R )Q]+Y J J dT - 
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Here fj are the Pauli matrices in Nambu-Gor'kov space. 
The trace operation implies summation over all possi- 
ble variables including replica indices and integration 
over energies. The equilibrium superconductive matri- 
ces Ql.r for the leads are diagonal in both energies 
and replicas and have the form Ql,b.( £ ) = t% cos^s^ ^ + 
T2 sin ipl,r in Nambu-Gor'kov space with <Pl(r) being 
the superconductive phase of the left (right) lead. This 
expression is valid at energies well below A. The contri- 
bution from higher energies |2j has already been taken 
into account by renormalization of the capacitance: 
C^C + e 2 {G L + G R )/(2A). 

To exclude fast fluctuations due to shifts of the elec- 
tron band by the potential <f> from the matrix Q we 
perform the following change of variables |S] 

0»(r) = K a (r), Qf T , = e *** a (T)ga^ e -w*JcV). 

(2) 

The phase K is determined up to a constant, which will 
be fixed later to simplify further analysis. With new 
variables the action takes the form 



S[Q,K] = --Tr(sr 3 Q) 



l/T 



dr 



4£ c 



2vE g (y) 



Q£> cos 2K + sin 2K 



(2), 



(3) 



In this formula we put ip = (pi, — ipu, denote Q^' — 
tr(fiQ)/2, N — C g Vg/e and introduce bare phase- 
dependent minigap 

5 
4 



E g (<p) = ^Gl + G 2 R + 2G L G R cos<p. (4) 



The expression (0 is very similar to the action for an 
SIN system with one superconductive lead . The only 
difference is the phase dependence of E g . Further cal- 
culation will essentially follow the procedure of Ref. PP . 
The key idea is the adiabatic approximation based on the 
separation of characteristic frequencies of matrix Q and 
phase K ensured by the inequality Ec S. Charac- 
teristic timescale of the variable K fluctuations is much 
shorter than that of electronic degrees of freedom, thus 
we integrate the action over K(t) regarding Q as a time- 
independent matrix (it depends only on the difference 
of its two time arguments). Than we apply the saddle 
point approximation to the if-averaged action. The jus- 
tification of this approximation will be provided below. 

We parametrize the time-independent matrix Q by 
an angle a: 



Q(e) — f 3 cos a(e) + fi sin a(e). 



(5) 



The T2-term here is eliminated by the proper choice of 
the constant in the definition of K. Inserting this ex- 



pression into @ we derive the Hamiltonian controlling 
the dynamics of the phase K: 

H = E C [{-id/dK - Nf - 2q{p) cos 2K] . (6) 

All physical quantities depend periodically on N. It is 
convenient to assume that \N\ < 1/2. The parameter 
q(tp) is expressed in terms of the angle a(e) 



I \ nEg(<p)T . 



(7) 



This sum diverges logarithmically and should be cut off 
at \e\ ~ A. For large values of q the phase K is nearly 
localized in the minima of cosine potential and the fluc- 
tuations are weak. In the opposite case fluctuations of 
the phase get strong and proximity effect is mostly sup- 
pressed. Thus the parameter q quantifies the strength of 
proximity coupling competing with the charging effect. 

With the derived Hamiltonian we are able to calcu- 
late the free energy of the K degree of freedom F(q, T) 
and then extract the total free energy of the system from 
the action © 

2nT 

T = —Y,£nCOsa(s n ) + F{q,T). (8) 

The equilibrium value of a(s) is determined by the mini- 
mum of this free energy functional (saddle-point approx- 
imation). This gives tana(e) = E g /e with E g obeying 
the self-consistency equation 



E a 



1 dF 
2Ec~dq' 



(9) 



This Eg is the minigap appearing in the spectrum of 
the normal grain. The estimation of matrix tQ fluctua- 
tions justifies the saddle-point approximation provided 
Eg » S. The system of two equations (J7J) and © de- 
termines two parameters, q and E g . A trivial solution 
q = Eg = always exists. It is analogous to the nor- 
mal state which is the stationary point (local maximum) 
of the superconductor free energy. We are looking for 
a non-trivial solution leading to nonzero value of the 
minigap E g . Once this solution exists it has lower en- 
ergy than the trivial solution. 

After solving the equations we can calculate the free 
energy © and all physical properties of the junction. 
We are interested in the current-phase relation given 
by the identity I(cp) = (2e/h)(dJ r /dip). Using the self- 
consistency relation @ and the identity Q we express 
the current as 



eS 2 E c E g q . 
= ^HEfW GLGRSmip - 



(10) 
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Below we consider analytically two limiting cases of 
weak ((ji> 1) and strong (9 « 1) charging effect and 
then discuss the numerical results for arbitrary q. The 
spectrum of the Hamiltonian © is given by the char- 
acteristic values of the Mathieu equation which is ele- 
mentary solved in these two limits. We first calculate 
the current at zero temperature taking the ground state 
of JBJ) for the free energy F. 

Weak Coulomb blockade. When charging effects 
are weak and the parameter q is large the phase K 
is localized near or 7r in the minima of the cosine 
potential. The applied gate voltage is ineffective in 
this case. Expanding the potential to the second or- 
der near its minimum we find the ground state energy 
Eo = Ec{—2q + 2^/q). The pair of equations (|7I9I) can 
be solved iteratively. In the considered regime the mini- 
gap is slightly suppressed in comparison with its bare 
value E g ((p). First, we estimate q substituting E g (ip) in 
the r.h.s. of 



<7o = 



E C S 



log(A/£ s0 )- 



(11) 



Here we neglect the (^-dependence of E g in the argu- 
ment of logarithm. At the next iteration we put q in 
the r.h.s. of © and then refine the value of q inserting 
the calculated E g into 0: 

E g =E g (<p)(l-l/2y/to) ) g = «fo(l-l/2V*). (12) 

With this q and E„ we calculate the Coulomb correction 
to the current using i|10|) 



I{<p) = I (<p) 



1 



1 



—G L G R log(A/E g0 )sin(p. 



(13) 



(14) 



Here we denote current in the absence of Coulomb in- 
teraction by Io(<p). Suppression of the current by the 
Coulomb interaction becomes stronger as the phase on 
the junction increases. Qualitatively, the bare minigap 
E g (ip) decreases and the charging effects win the com- 
petition with proximity further suppressing the current. 
In the symmetric junction (Gl = Gr) the proximity ef- 
fect vanishes completely as ip approaches n. The weak 
interaction approximation becomes invalid in this limit 
even if it is correct for small tp. 

Strong Coulomb blockade. For small values of the 
parameter q we calculate the ground state of 10 pertur- 
batively: E = -E c q 2 /2{1 - N 2 ). The equations lf7l5)l 
give 



Eg = 2 A exp 



q = 



2E g {l~N 2 ) ^ 



E g {<p) 



(15) 
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Fig. 1. Current vs. phase for the symmetric junction. 
Solid curve illustrates the case Ec = 0. It is given 
by the expression 1141 . Other curves correspond to 
Ec8/Eg = 0.5, 1.5, 2.5. At small values of this pa- 
rameter and at small tp the current is given by 1131 . 
At larger ip charging effects are always strong and the 
current is exponentially suppressed lib! . We assume 
Gl = G r = 20 and log(A/£ g0 ) = 5. 




J g0 



Fig. 2. The dependence of the critical current on 
Ec for the symmetric junction. The parameters are 
G l = Gr = 20, log(A/% ) = 5. 



Exponentially small minigap survives at T = when 
Coulomb blockade is strong. Josephson current is expo- 
nentially small as well: 



m = 2 4£^G L G R{ i 



KE${<p) 



N 



x exp 



4E C S 



(l-N 2 ) 



sin ip. (16) 



We solve numerically the system of equations I|7I9|I 
and plot the dependence of the current on the phase 
difference in Fig. Critical current I c = max^, I(ip) as 
function of charging energy is shown in Fig. [3 

The gate voltage enhances both the minigap and the 
current [see Fig- E] ■ Large Coulomb energy makes the 
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0.5 

C g V 3 /e 

Fig. 3. Relative increase of critical current by gate volt- 
age for Ecd/EgQ equal to 1.5, 2, and 2.5. This effect 
gets stronger at larger charging energies. Other parame- 
ters of the junction are Gl = Gr = 20, log(A/iJ 9 o) = 5. 




T5/E 2 g 

Fig. 4. Critical current vs. temperature for Ec5/Eg = 
0.5, 1.5, 2.5. Other parameters of the junction are 
G L = G R = 20, log(A/£ 90 ) = 5. 



charge of the grain to be nearly conserving quantity. 
Ground state corresponds to zero charge and is sep- 
arated by the gap AEc from the excited states with 
charge ±2e. States with odd charge are ineffective be- 
cause electrons tunnel from leads by pairs. This situ- 
ation changes when the gate voltage approaches e/2C. 
The gap between ground state and excited state gets 
twice smaller assisting tunneling of Cooper pairs. At 
higher gate voltage the ground state carries odd charge 
and the critical current starts to diminish. The increase 
of current with gate voltage is analogous to that studied 
in 0] where the similar setup with the superconductive 
grain was considered. 

Now we turn to the thermodynamic properties of 
the junction. The temperature dependence of the criti- 
cal current is found numerically and depicted in Fig. ^ 
At some temperature both the minigap and the Joseph- 
son current disappear. As temperature approaches its 
critical value T c the parameter q becomes arbitrarily 
small. This allows to expand the free energy of the K 
degree of freedom: F(q, T) = F{0,T) - E c (3{N,T)q 2 . 
The coefficient /3(N, T) may be found with the help of 
perturbation theory for the Hamiltonian pjl. 



(3(N,T) 



E- 



-(n-NfE a /T 



/(l-(n-A) 2 ) 



e -(n-N) 2 E c /T 

T « E c - 



2£ 

'l/[2(l-AT 2 )], 
E c /T - {2/3){E c /T) 2 , T » E c . 



(17) 



Note that (3 is non-monotonous function of tempera- 
ture. At large temperature highly excited levels of the 
Hamiltonian insensitive to the g-perturbation play the 
main role. Thus (3 goes to zero in this limit. At small 
temperature the phase K is almost frozen at the ground 
state. The g-term mixes two lowest excited states. As 



temperature grows these two states begin to contribute 
to the free energy increasing its dependence on q. When 
N approaches 1 /2 the ground state becomes degenerate 
and (3 falls monotonously with temperature. 
The expression (JJJ in the limit E g gives 



1 



EgEg 2 7 A 

log ■ 



E C S 



irT 



(18) 



Here log 7 ~ 0.577 is the Euler constant. The self- 
consistency condition in the limit of small q takes 
the form E g = E g (ip)(3q. Substituting this equation 
into Eq. I|18|l we find for the critical temperature: 



■2l 



A exp 



En 8 



Efa)0(N,T c 



(19) 



The same equation may be obtained by the expansion 
of (JSJ in powers of E g and setting the coefficient of E g 
to zero. It can be checked that the forth-order term of 
this expansion always remains positive. This justifies 
our assumption that E g vanishes continuously at the 
critical temperature. 

In the regime of strong Coulomb interaction critical 
temperature appears to be much less than Ec- Tak- 
ing low temperature asymptotic of [3(N, T) we find that 
T c = (j/TT)E g {T = 0), where E g (T = 0) is given by 
Eq. {TSJ. This is the BCS relation between the gap 
and the critical temperature. The phase of the grain 
strongly fluctuates and is mainly independent of the 
phase in the leads. The only effect of superconductive 
leads is a weak effective attraction in the Cooper channel 
which leads to formation of very weak BCS-like state. 

In the opposite limit of weak Coulomb interaction 
we employ high temperature asymptotic of (3 and find 



E 2 g ( v ) 2jAS 



log 



2 



(20) 
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Fig. 5. The critical line demonstrating the dependence 
of critical temperature on Ec- The three curves are 
plotted for AS/ Eg = 5, 17.7, 70. At large enough values 
of A this dependence may not be single- valued. Dot- 
ted line shows the turning points of critical tempera- 
ture. Small disk denotes the critical point where the 
dotted line touches the T c vs. Ec curve. This critical 
point gives the following critical values: T*5/E 2 = 0.54, 
E C S/E 2 = 2.05, A*5/E 2 g = 17.7. 

The whole dependence of the critical temperature 
on parameters is shown in Fig. It is possible that 
the equation (|19l) has more than one solution at a given 
value of Ec- This implies reentrant behavior of the 
minigap and critical current as functions of tempera- 
ture. Mathematically, it is due to nonmonotonic be- 
havior of the function (3{T). Physical explanation is 
most simple in charge rather than phase representa- 
tion of the Hamiltonian Two excited states with 
charge ±e have equal charging energies. Tunneling of 
Cooper pair, that switches between these two states, is 
not blocked by Coulomb interaction. Finite tempera- 
ture may excite the system to one of this states leading 
to the temperature-assisted proximity effect and the en- 
hancement of the minigap. The Hamiltonian JSJ con- 
serves the parity of electron number. Thus thermal- 
ization of the system implies some parity-breaking pro- 
cesses (e.g. single electron tunneling with energy above 
A), that may take a long time. 

At large enough values of A > A* = 17.7£^/(5, reen- 
trant behavior of the minigap with temperature was 
found in some (dependent on A/ A* ratio) interval of 
the Coulomb parameter EcS/Eg, cf. Fig. Fine tun- 
ing of the parameter EcS/Eg can be achieved by an 
appropriate phase bias, cf. Eq. 10} . 

To conclude, we have described the Josephson effect 
in a nanoscale SINIS junction modified by the Coulomb 
interaction. The most important feature of SINIS struc- 
ture (in comparison with usual SIS structure) is that 
it can demonstrate both good metallic conductance in 
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Fig. 6. The minigap in the normal dot as a function of 
temperature. Solid line is plotted for EcS/E 2 = 2.6. 
Dashed line illustrates the case E C S/E 2 = 2.67. Chain 
line is for E c 8/E 2 = 2.75. All curves are plotted for 
A = 70E 2 /S. The critical temperature as a function of 
Ec for the same value of A is shown by the chain line 
in Fig. U 

the normal state and Coulomb blockade of Josephson 
current at very low temperatures, since both the con- 
ditions Gl,r ^> 1 and Ej = J^I C < Ec can be ful- 
filled simultaneously. We calculated the current-phase 
characteristic of the junction in both weak and strong 
Coulomb blockade limit. The enhancement of the cur- 
rent by the gate voltage is predicted. The temperature 
dependence of the critical current and of the minigap 
induced in the normal part of the junction was found. 
A grain of noble metal with size about 50 nm connected 
to Nb superconductive electrodes by tunnel oxide barri- 
ers with transparency per channel of the order of 10~ 5 
could present an example of the studied system with 
Ej ~ E c ~ IK. We are grateful to T. Kontos and 
Ya. Fominov for useful discussions. This research was 
supported by the Program "Quantum Macrophysics" 
of the Russian Academy of Sciences, Russian Ministry 
of Science and RFBR under grant No. 04-02-16348. 
P. M. O. acknowledges financial support from the Dy- 
nasty Foundation and the Landau Scholarship (KFA 
Juelich). 
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